High-speed phonon imaging using frequency-multiplexed kinetic inductance detectors 
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We present a measurement of phonon propagation in a silicon wafer utilizing an array of frequency- 
multiplexed superconducting resonators coupled to a single transmission line. The electronic readout 
permits fully synchronous array sampling with a per-resonator bandwidth of 1.2 MHz, allowing sub- 
fis array imaging. This technological achievement is potentially vital in a variety of low-temperature 
applications, including single-photon counting, quantum-computing and dark-matter searches. 




Superconducting microwave resonators are a promis- 
ing detection technology for a variety of low-temperature 
applications. Electromagnetic resonators have been cou- 
pled to quantum information circuits^ and nanomechan- 
ical systems 2 -. They are also easily integrated with 
an absorber or an antenna, enabling them to be used 
for photon-limited astronomy, x-ray spectroscopy, or in 
dark-matter searches 3 " 1 - Their tunable resonance fre- 
quencies, narrow bandwidths and compatibility with 
large-bandwidth cryogenic amplifiers also admit the pos- 
sibility of frequency-domain multiplexing. This wire- 
sharing technique is widely recognized as being crucial 
to realizing large array sizes since the number of wires 
accessing the coldest stages of a cryostat is necessarily 
limited 5 . Frequency-domain multiplexing of a resonator 
design known as the Kinetic Inductance Detector (KID) 6 
has already been demonstrated at low read-out speeds 
and used in mm-wavelength astronomy^—. 

Despite their utility, the full potential of superconduct- 
ing resonators has remained unrealized. Single-resonator 
measurements typically use fast commercial electronics 
which give a large measurement bandwidth but are not 
compatible with array read out. On the other hand, cus- 
tomized array readouts utilizing programmable digital 
logic have not achieved per-resonator bandwidths exceed- 
ing ~100 Hz. Here we present the results from a recent 
experiment which realizes a large 1.2 MHz per-resonator 
bandwidth for an array of KIDs. The array is used to 
resolve the propagation of phonons created by the in- 
teraction of a cosmic-ray with the underlying substrate. 
Direct application of this phonon imaging sensor is im- 
mediately useful for cosmic-ray detection^ -, dark matter 
studies 1 ^ and for characterizing phonon propagation in 
exotic materials 12 . Further, by surmounting the techno- 
logical hurdle of array readout with a large per-resonator 
bandwidth, this measurement lays the ground work for 
developing, for example, a scalable quantum-computing 
architecture or single-photon camera based on supercon- 
ducting resonators. 

The operating principles of KID read out has been pre- 
viously discussed 6 . Briefly, the resonant frequency of an 
electromagnetic resonator is principally set by the geo- 
metric inductance and capacitance. In a superconductor 
however, energy can not only be stored in the electro- 



FIG. 1: (a) Array micrograph. Each line of five pixels was de- 
signed to occupy a 50 MHz bandwidth. The five labeled pixels 
were used to measure phonon propagation, (b) 2 mm x 2 mm 
LEKID. (c) In-phase (I) and Quadrature (Q) transmission 
amplitudes. The numeric labels correspond to the pixels in 
(a), the arrows indicate the direction of increasing frequency, 
and the star (*) marks the origin, (d) Total transmitted am- 
plitude for the measurement shown in (c), calculated using 
Amplitude 2 = I 2 + Q 2 . 



magnetic field, but also in the kinetic energy of the su- 
perconducting pairs. The resulting reactance is known as 
the kinetic inductance Lx ■ A superconducting resonator 
designed to maximize Lk will be sensitive to the inter- 
nal superconducting state. Typically, KIDs are designed 
to absorb photons or substrate phonons with energy ex- 
ceeding the superconducting gap. These can then break 
superconducting pairs and produce a measurable signal. 

For the current measurement, a KID geometry known 
as the the Lumped Element Kinetic Inductance Detec- 
tor (LEKID) was utilized 13 . A LEKID consists of a long 
inductive meander connected to an interdigitated capaci- 
tor used to tune the resonant frequency. An advantage of 
this geometry is that the current in the inductive section 
is uniform, resulting in a homogenous response to ab- 
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sorbed phonons throughout the meander. A micrograph 
of an array is shown in Fig. [TJa). To fabricate this array, 
a single 40 nm evaporation of aluminum deposited on 
270 /im thick silicon. The him was then patterned using 
standard UV lithography followed by wet etching. Af- 
ter dicing, the sample was affixed inside a superconduct- 
ing enclosure and connected to 50 fl coplanar-waveguide 
with wire bonds. The sample holder was subsequently 
mounted in a 100 mK dilution cryostat. 

The measured In-phase (I) and Quadrature (Q) trans- 
mission amplitudes for a 50 MHz frequency sweep of the 
array is shown in Fig. HJc) . The total transmitted am- 
plitude (IS21I) is shown in Fig. HJd). From this plot, the 
average loaded quality factor Ql is 1.3xl0 4 , limited by 
the coupling. By decreasing the coupling, can be in- 
creased up to the internal material-limited quality factor 
Qi (~ 10 5 ). This affords an increase in sensitivity at the 
expense of a decreased dynamic range. For the current 
measurement, a large dynamic-range was desired in or- 
der to differentiate between small and large events. The 
resonators were thus somewhat over-coupled. 
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FIG. 2: Electronic readout. A Digital-Signal Processing 
(DSP) card is used to frequency-multiplex the resonators. 
The individual components shown are: i) high-frequency syn- 
thesizer (.1-8 GHz), ii) splitter, iii) mixer, iv) attenuator, v) 
amplifier, and vi) low-pass filter. 

Programmable digital electronics were used to perform 
a fully synchronous transmission measurement of the ar- 
ray. The details of this technique have previously been 
discussed 8 . For the current measurement, a 12-bit, 100 
MSPS ADC was used for sampling the transmitted wave- 
form. This was followed by digital mixing, low-pass fil- 
tering and decimation resulting in a pair of I and Q val- 
ues for each resonator. Due to the slow roll-off of the 
digital 1.2 MHz low-pass filter implemented, an inter- 
resonator spacing of at least ~4 MHz was necessary to 
eliminate potential interference. Along with the 50 MHz 
Nyquist bandwidth of the ADC, the current electronics 
were thus able to measure up to 12 resonance in paral- 
lel. Using a faster ADC and improved digital filter, >100 
pixels per transmission line should be achievable with 
modern electronics. In order to reduce the high data- 
rate resulting from fast sampling, an individual trigger 
for each resonance was established. On-board memory 
simultaneously captured all of the signals if any trigger 
was exceeded. A diagram of the measurement setup can 
be seen in Fig. [2] 

The resonant frequency for an unloaded LEKID is 



fo 



2TT^(L K +L G )C(e) 



(1) 



where e is the effective surrounding permittivity, and 
Lk, Lq and C are the resonator's kinetic inductance, 
geometric inductance and capacitance respectively 1 ^. A 
change in the resonant frequency Sfo may be expressed 
as Sfo = SL K df /dL K . From this, Sfo = {-Cf^2)SL K . 
A small linear shift in the kinetic inductance results in 
a proportional frequency shift of the resonance feature. 
In the IQ plane, a small shift in frequency primarily re- 
sults in rotation around the resonance curve. This can 
be expressed as a rotation angle 4> as shown in Fig. [21a,) . 
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FIG. 3: Calibration, (a) The phase <j> as measured around the 
resonance curve in the IQ plane. The solid curve is a circular 
fit to the curvature near the resonance. The start (□) and end 
(o) points of the frequency sweep in (b) are displayed, (b) 
Frequency-dependence of <j> for a 120 kHz sweep around the 
resonance frequency. 

Since all the resonators are geometrically identical, ex- 
cept for slight differences in the capacitive section intro- 
duced to adjust the resonance frequencies, they all posses 
the same Lk- Hence, for a fixed energy input, all the res- 
onators will experience the same frequency shift A/. In 
Fig. [2Kb), a plot of <j> is shown for a 120 kHz frequency 
sweep around a resonance. A fit to the central data in 
this plot yields the slope A0/A/. For small cj), dividing 
by A<fi/Af results in a corresponding A/ and is effec- 
tively a calibration of the resonance. 

During a typical measurement, hundreds of events were 
captured over the course of a few hours at a rate of one 
event every ~5-20 seconds. An event causing a moderate 
response from all 5 pixels is shown in Fig.UJ The decreas- 
ing amplitudes in going from pixel 1 to 5 along with the 
increasing time delay corroborates the hypothesis that 
the cosmic rays are not interacting directly with the res- 
onators but rather create phonons in the substrate which 
propagate to the LEKIDs. At T=100 mK « T c -1.2 K, 
thermal phonons do not have sufficient energy to break 
superconducting pairs and cannot result in a measurable 
signal. The phonons which are detected are thus non- 
thermalized and posses energies exceeding the gap en- 
ergy. 

Collision-induced phonon propagation in silicon has 
been extensively studied . 14 ' 15 In the simplest model, nu- 
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FIG. 4: (a) Calibrated high-speed detection of a cosmic event 
using 5 KIDs. The individual responses are offset by 100 kHz. 
The dashed vertical line is the time that a resonator trigger 
level was exceeded, in this case by pixel I. (b) Zoom of initial 
rise clearly showing the propagation delay. 



clear or electronic recoil from a particle interaction cre- 
ates a cloud of high energy phonons (>5 THz). Rapid in- 
elastic scattering results in anharmonic down-conversion 
until reaching a bottleneck frequency of ~1.6 THz at 
which point they propagate outward, undergoing fre- 
quent elastic scattering. This quasi-diffusive processes is 
isotropic due to rapid conversion between phonon modes. 
The time difference of arrival between pixels and the fre- 
quency shifts thus enables a determination of the ini- 
tial event position and the propagation velocity. For the 
event in Fig. HI the position was determined to be at 
(-6.3±.3 mm,±3.6±.3 mm), where (0,0) is located at the 
geometric center of the array and the x-axis lies along the 
pixel centers. The degeneracy in the y-coordinate is due 
to the symmetry of the array which can be easily bro- 
ken in future measurements. The positional error is due 
to the arrival time uncertainty, equal to the inverse per- 
resonator bandwidth, and the ~2 kHz noise in the mea- 
sured frequency shift. Along with the thin substrate and 



complicating surface effects, the positional error on the 
order of the wafer thickness obscures any potential depth 
information for the current measurement. The propaga- 
tion velocity was determined to be 2.2±.9 km/s. 

By fitting the low frequency-shift tails of the cosmic 
events to an exponential, a time constant of 36 /is was 
measured. This significantly exceeds the resonator time 
constant r = (Ql)//o — 7 M s ) which is likely due to the 
long quasi-particle decay time in aluminum at T«T C . 
Finally, the initial rise of the event is caused by a para- 
metric shift in the resonant circuit which is not limited 
by r. The required resolution of the leading edge thus 
determines the minimum per-pixel bandwidth for a mea- 
surement. For example, in the current measurement, in 
order to resolve a ~ fj,s time delay between pixels, a per- 
pixel bandwidth exceeding 1 MHz was critical. 

In conclusion, we have imaged phonon propagation 
with sub-/is time resolution using frequency-multiplexed 
superconducting resonators. A measurement is planned 
in the near future utilizing a calibrated particle source 
in order to asses the detector sensitivity. An improved 
digital low-pass filter design, a larger bandwidth ADC, 
smaller pixels, and an increased pixel count will also be 
implemented in order to decrease the measurement un- 
certainty. Based upon these results, a full-scale detec- 
tor could then be developed for dark matter searches or 
cosmic-ray detection. This technique should also prove 
invaluable in developing a scalable quantum-computing 
architecture or single-photon camera based on supercon- 
ducting resonators. 
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